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Abstract 

In this paper the Casimir energy of two parallel plates made by materials of different penetration 
depth and no medium in between is derived. We study the Casimir force density and derive 
analytical constraints on the two penetration depths which are sufficient conditions to ensure 
repulsion. Compared to other methods our approach needs no specific model for dielectric or 
magnetic material properties and constitutes a complementary analysis. 
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I. INTRODUCTION 

One of the striking features of the Casimir effect jl, 2] is the dependence of the sign of the 
energy, which may be positive or negative, and of the corresponding repulsive or attractive 



staining 
4]. For 



force, on the geometry of the device and on the materials used. The possibility of o 
a repulsive force would open vast possibilities in the design of MEMS and NEMS 
example one of the principal causes of malfunctioning in MEMS is stiction: the collapse of 
nearby surfaces, resulting in their permanent adhesion. The possibility of having a repulsive 
Casimir force is an interesting way to avoid such a collapse of the structure but up to now 
there are only few experimental evidences for a repulsive force js-7] 

In this paper we study the Casimir energy of two parallel plates made by materials having 
different penetration depth. A study of the Casimir force and the role of surface plasmons 
between dissimilar mirrors has been carried out in the past for specific models of the dielectric 
and magnetic properties of the materials |8|. Here we propose an alternative method which is 
model independent and gives thus complementary information on the possibility of repulsive 
Casimir forces. We find that, depending on the relation between the penetration depths and 
the distance of the plates the force can be both attractive and repulsive. The penetration 
depth of materials can be taken into account by means of its connection to the surface 



impedance |9Nl lj . The surface impedance Z of any planar surface may be defined as the 
ratio of the complex electric and magnetic tangential field components at the surface 

E t (z ) = Z(H t {z ) x ft) (1) 

where n is a normal vector pointing inside the surface and Zq is the position of the surface. 
The main advantage of this formula is that it relates the tangential fields outside the material, 
thus it is not necessary to consider the internal degrees of freedom of the material which 
are taken into account through the values of Z [9]. Equation (prj) can be seen as an exact 
: unctional definition of the surface impedance so that it can be applied to arbitrary materials 



12] and it still holds when a description in terms of dielectric permittivity cannot be given 
131 ] . Indeed a complete correspondence with reflection coefficient and surface impedance 
exists 0, [3]. Moreover "/or large permeability and permittivity, the transition from 
attractive to repulsive behavior depends only on the impedance Z" 

The paper is organized as follows: in section II the Casimir force in the general configu- 
ration is evaluated and some limiting results are recovered. In section III the conditions for 



having repulsions are derived. Finally section IV contains remarks and conclusions. 



II. THE CASIMIR FORCE 

In the following we will consider two parallel plates lying in the (x, y) plane, located at 
z = 0, z = a and characterized by different surface impedances Z(p t a)(u>) respectively. Given 



the functions <5(o, a )(w 



•• Z (0,o)(w) 



material at the frequency lo 



10 



Re[8{uj)] is interpreted as the penetration depth of the 



llj see also 13( . Because of translational invariance in the 



(x, y) plane the electric and magnetic fields can be written as (in the following we will use 
natural units: h = c = 1): 



E(x, t) = f{z)e ik ^ x± - ut , B(x, t) = ^{z)e ik ^- ut 
with k± = (k x , k y ) and x± = (x,y). The Maxwell equations imply: 

0: 



_ + A-f = 0;J + A-, 



(2) 



(3) 



with A 2 = uj 2 — k 2 x — k 2 . Imposing relation ([TJ) we obtain, for the (x, y) components of E and 
B the following boundary conditions [ill Il3| . 



(4) 



/.(0) = -5 (ik x f z (0) - f>(0)) f y (0) = -5 (ik y f z (0) - f y (0)) 

fx(a) = $a(ik x f z (a) - f x (a)) [ f y (a) = S a (ik y f z (a) - f y (a)) 

moreover everywhere V • E = must be satisfied. 

In this way, with a suitable choice of the reference frame, we find the following dispersion 
equation: 

A TM = [5 5 a (X 4 + k]_) + X 2 (25 5 a kl - 1)] sin (aX) - (5 + 5 a )X(X 2 + k 2 ± ) cos (aA) = (5) 

for the TM modes and 

A TE = 5 5 a (X 2 - 1) sin (aA) - (<J + 5 a )X cos (aA) = (6) 

for the TE ones. We use the argument theorem to obtain the Casimir energy so that, 

after u- rotation to the ima gina ry axis: u> — > i(, we have (for the properties of S (or Z) along 
the imaginary axis see 10, llll, ll5|] ) 

1 



E 



d(dk x dk y In [A TM (A, i()A TE (X, 



(7) 



This integral diverges and, as usual, to regularize it we must subtract the energy correspond- 
ing to the configuration with the two plates infinitely far away (a — > oo): 



A°c • an ( l X X 2 5 0<? ^ a q V 



A OO 



ie 



aq 



(d q 2 -q- 5 k x 2 ) (5 a q 2 -q- 5 a k x 2 ) 



with q = \J( 2 + k\. Thus the renormalized Casimir energy will be given by: 



2(2tt) 



ln[ 1 - e 



d(dk x dk v In 



(l + ^ a (iC))(l + ^0«))/ 



+ 



-2aq 



(g-^«)C 2 ) (q-So^OQ Y 

(q + ^C)C 2 ) (g + 5o«)C 2 )J J 



or, in dimensionless variables 
1 



E 



R 



oo ;>oo 



d£ / k±dk± In [ 



(l+p5 a ) (l+p6 ) 



+ 



Air 2 a 3 Jo 
ln[( l-e~ 2p 



with p = ag, fc (X)W) = ak { ^ y) ^ = a(, 5 {afi) = 

In the following we will concentrate on the Casimir force, it can be written: 



(p - 5„£ 2 ) 


(p - s e) 


(P + 8 a ?) 


(p + U 2 ) 



1 



2ir 2 a 4 Jo 



OO POO 



d£ p dp 



- 2 p(5oP-1)(Lp-1) 



(5 oP + l)(5 a p +!)(!- 



(<5 OP +l)(5aP+l) 



+ 



(8) 
(9) 



(10) 



(11) 



(12) 



Now it is not difficult to show that the contribution coming from the point £ = is zero, thus 
we can safely remove this point from the integral, which allows us to rewrite the integral: 



'R 



| roo roo 00 



2pn 



27r 2 a 4 

— 2pn 

00 
n=l 



(6 oP -l)(6 aP -l) 
(5 p+l)(5 a p+l) 



+ 



{p - s e) (p - ie) 
ie+p)_ 



(13) 
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The absolute values of the terms in the brackets are always less or equal to one, in the case 
= S a = (0, oo) they are maxima (1) and we have: 



X I POD POD 

^ = -E^/ P^pe- 



2pn 



n=l 

oo 

-E 



E 

n=l 



7T 2 a 4 Jo 



e" 2n «(2<«+ 1) + 1) 
4n 3 



7T 



^ 8a W 



240a 4 



(14) 



In contrast, if we take 5q = oo, 5 a = or viceversa they take on minimum values of (—1) 
and we obtain 



°° I POD POO 



n — 2pn 



oo o 

E° 



7 7T 2 



a^n^vr 



8 240a 4 



16 



(15) 



for two non dispersive 



17] . The upper calculation 



Thus, in this case we recover the result obtained by Boyer 
mirrors having e = (oo, l),fi = (l,oo) respectively, see also jgj 
also constitutes an independent demonstration of the result found by Henkel and Joulain 
eq.(4) of 

From eq. ( fl3l) we may understand intuitively what kind of conditions must be satisfied 
to have repulsion. Indeed, if the two slabs are made of the same material we have 5 a = 5q 
and the expression of the force becomes: 



E 



\ 2vr 2 a 4 Jo 



d£ / p 2 dpe 



—2pn 



's p - r 

JoP+l, 



2» 



+ 



M 2 +p, 



2/i." 



(16) 



In this case the integrand is always positive and the force will be always attractive. The only 

(p-soe)( P ~Saey n 

(p+<5o£ 2 )(p+<U 2 ). 



(6pp-l)(S aP -l) 

(5 p+l)(5aP+l) 



+ 



possibility to have repulsion is to have S a ^ So, such that 
be negative. Fortunately, the series starts with the term n = 1 so that the possibility is not 
ruled out. 

In the next section we will study the case 5 a 1 and we will concentrate on the first 
term of the series: n — 1. 



III. SUFFICIENT CONDITIONS FOR REPULSION 



In the following we will develop F% at first order around 8 a — 0. After the integration on 
the p variable we will study the behavior of the remaining integrand which will be a function 



of £ only and determine what conditions must be satisfied to have repulsion. 

2S P \ 





/ 


2vr 2 a 4 


Jo 




1 


8vr 2 a 4 


Jo 


f 





d£ / p 2 dpe- 2p 



(2pS a - 1) 1 



SoP + l, 



p A ^e + WJ 



128vr 2 a 4 Jo 



4 3 . 

i 



(17) 



With: 



A(0 = -4£ 4 5 5 + (4£ 3 + 2? - 24 (2£(£ + 1) (2£ 2 + 3) + 3)) 5 4 + 

4(2£(4£(£ + 2) + 7) + 7)S 3 + (5 a (-8£ - 4) - 4£ - 2)5 2 + (84 + 4)4 

12(0 = 6^(0 



E n (x) is the exponential integral function 19) . Let us study /#(0 for the two regimes, 
< 4 < So <C 1 and < 4 -C 1, 4 ^ 1- In the first case we find: 
1 



3vr 2 a 4 



[(So + 6a) (8£ 3 + 8£ 2 + 6£ + 3) - (2f 2 + 2£ + 1 



(18) 



Thus, in the range of frequencies relevant to the Casimir effect, £ ~ 1, the condition for 
having repulsion is, at second order in £: 

2 + 401 + 



S > 



S n 



-S a + 0.775 - 0.494£ + 0.119O- 



(19) 



3 + 2^(3 + 4^(1 + 0) 

This shows that it would be possible to have repulsion if S a ~ and Sq > 0.4. However this 
last condition is in contradiction with the assumption 5 1. Let us also note that if we 
assume S a = S = const we can evaluate all terms of the series ( TTBl and recover the result 



of Mostepanenko and Trunov [11], (see also 18] for equivalent results for a scalar field). 
If < 5 a <C 1, S ^> 1 we find at first order (for the asymptotic expansion of E n (£) see 

F1r = 16+W / d ^ {-% A &>® + *> 5 &>0 + C &>fl} 

with: 



(20) 



A(5 a , = -240(3 + 2£(3 + £(3 + 2f))) 

5(4, = 2£ 2 + 40(1 + + 40(7 + 2^(7 + 60) 

C(S a , = -3 - 6£ - 80 - 80 + 840 



The term within the curly brackets in Eq. ( 1201) is a second order polynomial in So and, to 
have repulsion, it must be positive. Since the coefficient of the 5$ term is always negative 
we must require that the discriminant of the associated second order equation is positive. 
This discriminant is a second order polynomial in S a and we have to study the associated 
second order equation: 

%D(t)+5 a E{£) + F({) = (21) 



with 



£>(£) = 49£ 4 + 196£ 5 + 556£ 6 + 720£ 7 + 528£ 8 + 256£ 9 

E(£) = -72£ 2 - 288£ 3 - 596£ 4 - 848£ 5 - 744£ 6 - 432£ 7 - 160^ 

F(0 = 4£ 4 + 16£ 5 + 32£ 6 + 32£ 7 + 16£ 8 . 



Since D(£) > 0, for the polynomial to be positive we have to choose 5 a such that it lies 
outside the interval defined by the two roots 5 a nm and 5 a nax of the corresponding equation, 
that is l a < 5™ n or 5 a > l™ ax for G = E 2 -4DF > 0. If G < we may take any value for 8 a . 
In that case the two roots coincide if 5 a = and become imaginary of 5 a > excluding any 
physical solution. Fig.l illustrates the behavior of the two roots as a function of imaginary 
frequency 
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FIG. 1. The two roots of eq. (|2ip . 5 a nm and <5™ ax , are shown as function of imaginary frequency £ 
(solid and dashed curve respectively). For £ > 2.85, G < 0, at £ = 2.85 we have <5™ n = $<£ ax and 
at larger frequency any value of b a will give rise to a positive value for (|2ip . 
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Since 5™ ax can be larger than 1 (see Fig.l) and we assumed 5 a C 1 we remain with the 
only possibility S a < 5™ in '. Around £ ~ 1 the condition for the positivity can be written: 



S a < £™« = — - — ~ -0.0135 + 0.05148£ - 0.00533f 2 . (22) 

If this inequality is satisfied the force density will be repulsive for those values of 5 which 
satisfy 

$o <o <o , (23) 
fimin ^ fimax ^ \y e [ n g fa e sma ll er (larger) roots of the associated equation: 

- 6 2 A(5 a , + 5oB(L + C(5a,Q = (24) 

~ B + VB " ~ AAC a +lin 
<J = ~ ai + d a a 2 

o Q = — ~ a 3 - -^a 4 - d a a 2 

where 

ai = 10.48 - 12.56£ + 4.58£ 2 , a 2 = 152.438 - 230.856£ + 93.118£ 2 , 

a 3 = -9.16122 + 12.0305£ - 4.50083£ 2 , a 4 = 0.376294 - 0.126549£ + 0.013413£ 2 

Note that in the case of an ideal mirror at z = a we have 5 a = 0, 5™ ax — > oo and the only 
condition to be satisfied is: 

5 > ~5™ in = 10.48 - 12.56£ + 4.58£ 2 . 

Thus the situation in which one mirror is ideal gives rise to quite different results than 
the ones obtained when both are real. When both mirrors are real, they both must satisfy 
restrictions to ensure repulsion and, moreover, a precise relation between the two penetration 
depths must be fulfilled ( eqs. (1221231) ). 

IV. DISCUSSION AND CONCLUSIONS 

Let us apply our results to one known situation, that is of two mirrors described by the 
plasma model. In this case we have 

$(o a) = ; with = auf n n 

/T^p Y _|_ £2 s (0,a) (0,a) 



uf \ being the plasma frequency of the mirror in z = (0, a) respectively. Condition ( 122 



gives > ~ 30.65 which means that the mirror in a must have a plasma frequency 

up > 30.65/a, but condition (EH} implies , 1 > 4.37 which, being J{$) 2 + 1 > 1, is 

V W + 1 

impossible. 

Let us consider now the case of hypothetical materials having 8q(£) = with k = 
1,2,3,4. There we obtain for the Casimir force, using the exact first three terms of the 
series eq. (|T3|) : 

a 4 F R (k = 1) = -0.0007 - 0.0004 - 0.0000 = -0.0011 
a 4 F R (k = 2) = 0.0084 - 0.0005 + 0.0000 = 0.0079 
a 4 F R (k = 3) = 0.0133 - 0.0006 + 0.0001 = 0.0128 
a 4 F R (k = 4) = 0.0164 - 0.0008 + 0.0001 = 0.0157 

The result is illustrated on Fig. 2. The left hand part shows 8™ m (dashed line) and 5 a = 
y/30. 65 2 g 2 ( so ^ nne ) as a function of imaginary frequency. In the right hand part the 
shaded area gives the range of values of 5q given by condition ( 123]) for which the Casimir 
force becomes repulsive while the dashed, dotted, dotted-dashed and solid lines, decreasing 
monotonously with increasing frequency, represent 5q for k = 1,2,3,4 respectively. For 
k = 2,3,4 the force turns out to be the more repulsive the higher the fc-value, even though 
the values of 5q are only on the limit of the favorable region. 
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FIG. 2. In a 5™ n and 5 a 



are shown as dashed and solid curves respectively. In b the 



area between the two curves Sq 1 ^ , <5™ n (shaded area) and the curves for k = 1, 2, 3, 4 dashed, 
dotted, dash-dotted, and thick respectively are shown. For £ < 1 5 a > 5™ m and consequently the 
results are imaginary and no physical solution exists. 
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In conclusion we have found sufficient conditions for the Casimir force to be repulsive 
with an approach considering only the skin depth and needing no specific model of dielectric 
or magnetic properties. It would be interesting to study now how much these conditions 
can be softened, as after all to have a positive integral it is not necessary to have a positive 
integrand. From this point of view our analysis must be deepened trying to obtain analytical 
necessary conditions for having a repulsive force. Nonetheless our result demonstrate that 
repulsion is possible if the penetration depth of the two mirrors satisfy appropriate relations. 

The approach seems promising as it can be extended to anisotropic material characterized 



121. 



by a tensorial surface impedance and to more general material 

It would also be very interesting to derive the skin depth at optical frequencies from the 
available tabulated data to search for materials matching the conditions we have established 
and to use our result to design new materials such as to have repulsive properties. 
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